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ABSTRACT 

With the goal of investigating the nature of O vi absorbers at high redshifts, we study 
the effects of proximity to the background quasar. In a sample of sixteen quasars at 
z qso between 2.14 and 2.87 observed at high signal-to-noise and 6.6 kms -1 resolution 
. with VLT/UVES, we detect 35 O vi absorption-line systems (comprising over 100 indi- 

£L> vidual O vi components) lying within 8 000 kms -1 of z qso . We present component fits 

to the Ovi absorption and the accompanying Hi, Civ, and Nv. The systems can be 
£5 ■ categorized into 9 strong and 26 weak O vi absorbers. The strong (intrinsic) absorbers 

' are defined by the presence of either broad, fully saturated O vi absorption or partial 

coverage of the continuum source, and in practice all have logiV(0 vi)>15.0; these 
systems are interpreted as representing either QSO-driven outflows or gas close the 
central engine of the AGN. The weak (also known as narrow) systems show no partial 
coverage or saturation, and are characterized by log N(0 Vl)<14.5 and a median total 
qq ■ velocity width of only 42 kms" 1 . The incidence d7V/dz of weak O vi systems within 

2 000 kms -1 of the quasar down to a limiting equivalent width of 8 mA is 42±12. 
■ Between 2 000 and 8 000 kms' 1 , dN/dz falls to 14±4, equal to the incidence of in- 

' tervening O VI absorbers measured in the same spectra. Whereas the accompanying 

\^ , Hi and Civ column densities are significantly lower (by a mean of ~1 dex) in the 

f*^ ■ weak O vi absorbers within 2 000 kms -1 of z qso than in those at larger velocities, the 

OO | O vi column densities display no dependence on proximity. Furthermore, significant 

offsets between the H I and O vi centroids in «50 per cent of the weak absorbers im- 
ply that (at least in these cases) the H I and O vi lines are not formed in the same 
phase of gas, preventing us from making reliable metallicity determinations, and rul- 
ing out single-phase photoionization- model solutions. In summary, we find no firm 
evidence that quasar radiation influences the weak Ovi absorbers, suggesting they 
are collisionally ionized rather than photoionized, possibly in the multi-phase halos 
of foreground galaxies. Non-equilibrium collisional ionization models are needed to 
explain the low temperatures in the absorbing gas, which are implied by narrow line 
widths (b< 14 kms -1 ) in over half of the observed O vi components. 

Key words: cosmology: observations - quasars: absorption lines - intergalactic 
medium. 
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1 INTRODUCTION 

Quasars are ideal backlights for absorption-line spec- 
troscopy, thanks to their high luminosities and flat con- 
tinua. Yet in addition to allowing us to detect foreground 
absorbers, quasars also create and influence them, through 
the effects of outflows and ionizing radiation. The enhanced 



* Based on observations taken under ESO Programme ID 166. A- 
0106(A), using the Ultraviolet and Visual Echelle Spectrograph 
(UVES) on the Very Large Telescope, Unit 2 at Paranal, Chile, 
f E-mail: afox@eso.org 



level of ionizing radiation in the vicinity of quasars gives rise 
to the line-of-sight proximity effect, the observed decrease in 
the number density of Lya forest lines (equivalent to a de- 
crease in the mean optical depth, and an increase in the level 
of hydrogen ion ization) at velocities ap proachi ng the quasar 
( Carswell et all Il982l; iMurdoch et~afl Il986l: iTvtlerl Il987l: 
Baitlik Duncan, fc Ostrikerl 1 19881: iLu. Wolfe, fc Turnshekl 



199ll: iBechtoldl Il994l: iRauchl Il998l: IScott et all l200d : 
Guimaraes et al.ll2007l ; iFaucher-Giguere et al.ll2008h 



To obtain an unbiased view of absorbing structures 
in the high-z Universe, one needs to remove the effects 
of proximity from samples of metal-line systems. The 
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broad absorption line (BAL) systems, easy to identify in 
QSO spectra because of their optically thick absorption 
extending over thousands of kms -1 , are clearly connected 
to the quasar. BALs trace high-velocity QSO o utflows at 
speeds of up to 0.1c (e.g. iTurnshekl Il984 Il988l). and are 
seen i n 10-15% of quasar spectra at z > 1.5 (|Wevmann et al.l 
ll99ll ; lTrump et al.ll2006h . A separate category of mini-BAL 
systems, with fully saturated absorption profiles and total 
velocity w i dth < 2 000 kms -1 , has a l so been identified 
dTurnshekl 1 19881; IChurchill et al.1 Il999l ; lYuan et al.l |2002| ; 
iMisawa et all l2007bl ). " As with BALs, mini-BALs are 
interpreted as tracing outflows ejected by QSOs. A more 
difficult task is to decide whether narrow absorbers in 
quasar spectra are either created by or ionized by the 
QSO. The traditional way to do this employs the displace- 
ment velocity from the quasar, Sv = v qso — u a bs- So-called 
'associated' or 'z a bs ~ z qS o' or 'proximate' absorbers, 
typically defined as na rr ow absorbe r s at Sv < 5 000 kms -1 
JWeymann et al.l fl979l; iFoltz et al.l Il98rj ; lAnderson et all 



2000| ; ICarswell. Schave. fc Kiml |2002|; IBergeron et al 
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Hamann fc Ferlandl ll999F ~ are often removed from 



intervening samples because of the possibility they may 
be intrinsic^. Arguments in favour of an intrinsic nature 
for many proximate absorbers include time-variability 
of absorption, partial coverage of the continuum source, 
detection of excited ionic states implying high gas densities, 



suoer-solar metallicities (IWamDler. Bergeron. & Petitiean 


1993 


; Mailer. Jakobsen, & Perrvman 


1994 


; Petitiean. Rauch, & Carswelll 


1994; 


Savaglio, D'Odorico, & Moller 1994 iTriDD et al.l 


1996; 


Barlow & Sargentl 19971; Hamann 1997; Hamann et al. 


1997313. l200ll; Hamann. Netzer, & Shields 


2000; 


Petitiean & Srianandl Il999l; Srianand & Petitiean! 


2000; 


Ganeulv et al. 1999, 200ll. 120061; Naravanan et al.l 2004h. 



Using a velocity cutoff to differentiate inter- 
vening and intrinsic absorbers has two main prob- 



lems: quasar-cjcctcd intrinsic systems can 
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higher Sv IIHamann. Barlow, fc Junkkarinen 



app ear at 
1997bl; 



Richards et~al1 Il999l; iRichardsl l200ll; IMisawa et ail 



2007al : 

and intervening 



Nestor. Hamann. fc Rodr iguez HidalgoJ 2008h . 



system s can appear at lo wer Sv 



|M0ller. Warren, fc Fvnbol Il998l ; ISembach et ail |2004|) . 
So are many genuine intervening systems lost when they 
are excluded from absorber samples just because of their 
proximity to the quasar? In this paper we investigate the 
transition between proximate and intervening absorbers. 
We focus on O VI absorbers, which trace either warm-hot 
(T> 10 J K) collisionally-ionized plasma, or photoionized gas 
subject to a hard ionizing spectrum extending to energies 
above 113.9 eV (the ionization potential to create +s ). 
O VI absorbers are of interest for many reasons, including 
the significant role they play in the baryon and metal 
budgets, the window they provide on intergalactic metal 
enrichment, and their ability to trace energetic galaxy/IGM 
interactions such as accretion and galactic wi nds. After the 
first d etection of Ovi absorbers at z=2-3 bv lLu fc Savage! 
(1993), much work has been done in the era of lOm-class 
telescopes to characterize their properties l|Schave et al.l 



We favour the use of the term 'proximate' since it involves no 
assumption about the absorber origin. 



ISimcoe. Sargent, fc RaricrT 2002 7 120041 



Simcoe et al 
20071 ; 



2002 

1 20061; IBergeron fc Herbert- Fortl 120051; lLopez et al.1 
lAguirre et al.ll200Sl ; iGoncalves. Steidel. fc PettinilbOOcf T 

Here we ask a simple question: is there a proximity 
effect in O VI? In other words, can we see the signature 
of quasar photoionization, in terms of correlations between 
the ionization properties of absorbers and their proximity 
to the quasar? To address this issue, we form a sample of 
proximate O VI absorption systems at z=2-3 using a homo- 
geneous set of high-resolution VLT/UVES quasar spectra, 
and we then compare their properties to a sample of inter- 
vening (i.e., non-proximate) Ovi absorbers observed in the 
same set of spectra; the intervening samp le has been pub- 
lished by IBergeron fc Herbert-Fort! (|2005l . hereafter BH05), 
but it has since been enlarged from ten to twelve sight lines. 
We also briefly compare our results on Ovi absorbers at 
2=2-3 to those obtained in the low-redshift Universe with 
space-based ultraviolet spectrographs. This paper is thus a 
study of both the proximity effect and the nature of Ovi 
absorbers in general, and is structured as follows. In §2 we 
describe the data acquisition and reduction, and our ab- 
sorber identification and measurement processes. In §3 we 
discuss the observational properties of proximate Ovi ab- 
sorbers. In §4 we discuss the implications of our results, and 
we then present a summary in §5. Throughout this paper we 
adopt a W MAP 3-year cosmolo gy with f2M=0.27, f2A=0.73, 
and /i7o=l (|Spergel et al.ll2007l ). 



2 DATA ACQUISITION AND HANDLING 
2.1 Observations 

The ESO-VLT Large Programme The Cosmic Evolution 
of the Intergalactic Medium (IGM) has built a homoge- 
neous sample of high resolution, high signal-to-noise quasar 
spectra for use in studying the high-redshift IGM. Twenty 
bright quasars (most with V < 17) in the redshift range 
2.1-3.8 (median z qsc> =2.44) were selected for observation, 
deliberately avoiding sight lines that contain high-column 
density damped Lyct absorbers. The observations were car- 
ried out with the Ultraviolet an d Visual Echelle Spectro- 
graph (UVES; [Dekker et al.ll2OO0l ) on VLT UT2 at Paranal, 
Chile, using a total of thirty nights of service-mode tele- 
scope time between June 2001 and September 2002. The 
observations were taken in good seeing conditions (<0.8") 
and at airmass<1.4, with both blue and red dichroics (Dicl 
346+580, Dicl 390+564, and Dic2 437+860 settings), giv- 
ing spectra complete from 3 000 to 9 000 A. The data reduc- 
ti on was conducted usi ng the MIDAS pipeline described 
in lBallester et aL (2000), and then continua were fitted us- 
ing a third-degree spline function interpolated between re- 
gions free from abs orption. Full details of t he reduction 
process are given m IChand et all (|2004l ) and lAracil et al.l 
(2004). The 2x2 binning mode was used, giving a rebinned 
pixel size of 2.0-2.4 kms -1 . The spectral resolution of the 
data, R ~ 45 000, corresponds to a velocity resolution of 
«6.6 kms -1 (FWHM). Typical signal-to -noise ratios per 
pixel in the spectra are «35 at 3 500 A and «60 at 6 000 A. 
Among the twenty quasars, sixteen are at z < 3. These form 
the basic sample for this paper. The four cases at z > 3 were 
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excluded, because the density of the Lya forest at these red- 
shifts renders futile searches for Ovi. 

The emission-line redshifts of ten of the sixteen 
quasars in our samp le have been carefully determined by 
iRollinde et al.l (|2005l . and references therein) , and we adopt 
their values. For the remainin g six quasars we take z qS o 
from IScannapieco et all l)2006h . Each of these quasar red- 
shifts is based on a combination of the observed Ha, Mgn, 
Civ, and Siiv emission-line redshifts, allowing for system- 
atic shifts in the displac ement velocities of the various lines 
fe.g. iTvtler fe Fanlll992f). The uncertainty on each z qso is of 
the order 500 kms 1 l|Rollinde et al.ll2005h . 



2.2 Identification of Ovi systems 

We searched systematically through the sixteen quasar spec- 
tra in our sample for Ovi absorbers within 8 000 kms -1 of 
the quasar redshift. This velocity range corresponds to red- 
shifts between z m i n =0q SO -[(l+2qso) 8000 kms -1 /c] and z qso . 
Our choice of velocities up to 8 000 kms -1 from the quasar 
was made in order to investigate whether a transition occurs 
at 5 000 kms" 1 , the conventional cutoff between proximate 
and intervening systems. In practise, we also searched for 
(and found) Ovi absorbers at up to 2 500 kms" 1 beyond 
the quasar redshift, indicating the presence of high- velocity 
inflows toward the quasar. 

32 O VI absorbers were identified by their presence 
in both lines of the Ovi AA1031.926, 1037.617 dou- 
blet in the correct ratio through the line profile, i.e. 
T a (v, A1031)/r a (ii, A1037) = 2, where the apparent optical 
depth T a (v) — In [F (v)/F(v)], and where F(v) and F c (v) 
are the observed flux level and estimated continuum flux 
level, respectively, as a function of velocity. In addition, three 
Ovi absorbers were identified in one line of the Ovi dou- 
blet only, with partial blending in the other O VI line, but 
with corresponding absorption seen in C IV and H I at ex- 
actly the same redshift. Together the sample consists of the 
35 Ovi systems listed in Table 1. We define the system red- 
shift z a bs by the position of the strongest O VI component in 
the system, i.e. the wavelength of maximum optical depth. 
In each system, we looked for accompanying absorption in 
Hi (Lya, Ly/3, and Lyy), Cm A977, Civ AA1548, 1550, and 
Nv AA1238, 1242 at the same redshift z abs . 

We denote the velocity displacement of each absorber 
from the quasar as Sv = c|z qso — z a bs|/(l + z a bs). A nega- 
tive value corresponds to an absorber moving toward the 
quasar. The O VI systems almost always contain multiple 
components clustered together. We thus define Sv both for 
each absorber and for each individual Ovi component. The 
uncertainties on all velocity displacements are dominated 
by the ~500 kms" 1 uncertainty of the quasar emission red- 
shifts. 



2.3 Classification of Ovi systems 

After identifying 35 proximate Ovi systems, we found that 
they could be classified into: 

(a) Seven strong O VI systems, which are defined by 
showing both fully saturated O VI lines extending over tens 
to hundreds of kms -1 , and strong (often saturated) Nv and 
C IV absorption. In practice we find that these seven systems 



all show log JV(0 vi)> 15.0 and are at Sv < 3 000 kms -1 , with 
four of the seven at Sv < (z a b s >2; qso ), i.e. the strong sys- 
tems are strongly clustered around z qso . The strong systems 
are interpreted as being intrinsic to the quasar. Three of 
the strong systems are classified as mini- balU but no full 
(classical) BAL systems are present in our sample. Partial 
coverage of the continuum source is present in many of the 
strong systems, as evidenced by the failure to find a suc- 
cessful Voigt profile fit to the Ovi profiles, even though the 
velocity structure is common among the two members of 
the doublet. This confirms an intrinsic origin, because par- 
tial coverage of the continuum source only occurs when the 
absorbing gas is physically close to the AGN. The strong 
systems contain multiple components, with a mean of 5.9 
components per system. 

(b) 26 weak O VI systems, also known as narrow sys- 
tems, which are unsaturated and show no evidence for par- 
tial coverage, i.e. we were able to conduct successful Voigt 
profile fits to the two lines of the O VI doublet. The weak sys- 
tems all have 12.80<log N(0 vi)<14.50, clearly separated 
from the strong systems, but coincident with the range 
of O VI column densities observed in intervening systems 
(BH05). Their total line widths are narrow, with a median 
and standard deviation of A«go=42±33 kms -1 , where Avgo 
is the line width encompassing the central 90% of the inte- 
grated optical depth. The weak O VI systems are accompa- 
nied by Civ absorption in 24/26 cases, and in 9/26 cases by 
Nv. An mean of 1.9 components per system is required to 
fit these weak absorbers. 

There were two intermediate cases that presented diffi- 
culties in their classification: the systems at 2=2.6998 and 
2.7134 toward HE 0151-4326. Neither of these systems shows 
a saturated O VI profile, and their O VI column densities lie 
in-between the values shown by the strong (>15.0) and weak 
(<14.5 populations). However, differences between the pro- 
files of the two doublet lines of O VI, and between the two 
lines of N v, in similar velocity ranges are strongly suggestive 
of partial coverage. This implies they should be classified as 
intrinsic. For this reason we place these systems in the strong 
category, and our final sample consists of 9 strong systems 
and 26 weak systems. 

In addition to the different O VI (and N v and C iv) col- 
umn densities, there are several lines of evidence that sup- 
port the notion that the strong and weak samples trace phys- 
ically different populations. The strong population have ve- 
locities (relative to z qso ) that distribute differently, with 4/9 
at z a bs>-z qS o, whereas 22/26 weak absorbers are at z ab s ^ 2: q so - 
The strong population are show markedly different Hl/O VI 
ratios than the weak population (see §3.5). Finally absorbers 
in the strong population frequently show evidence for partial 
coverage of the continuum source. Whereas the strong ab- 
sorbers are clearly associated with the quasar (e.g . lTurnshekl 
1 1984 IWevmann et all Il99ll : lHamann fc Ferlandl Il999h . ei- 
ther tracing QSO outflows or inflows near the central en- 
gine of the AGN, the origin of the weak O VI absorbers is 
less clear. They are often called 'associated' because of their 
proximity in velocity to the quasar, but there is not neces- 



2 The mini-BALs are the absorbers at z a j, s =2.4426 toward 
HE 1158-1843, 2.1169 toward HE 1341-1020, and 2 abs =2.9041 
toward HE2347-4342. 
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sarily any physical connection. In the rest of this paper we 
investigate the properties of this weak, proximate sample. 



2.4 Component fitting 

We used the VPFIT software package^ to determine the col- 
umn densities, line widths, and central velocities of the O VI 
components in each system. The number of components to 
be fit was usually self-evident from inspection of the line pro- 
files, but in case of doubt we used the minimum number of 
components necessary. Our sample of 35 O VI systems com- 
prises 101 components. We also fit, when possible, the Hi, 
Nv, Cm, and Civ absorption, independently of the Ovi 
fits, i.e. we did not fix any component velocity centroids or 
line widths when fitting the various ions. In the case of the 
H I lines, we simultaneously fit Lya with whichever of Ly/3, 
Ly7, and Ly<5 were unblended. The total column densities in 
each ion in each system were obtained by summing the com- 
ponent column densities, with the individual column density 
errors added in quadrature to produce the overall error es- 
timate. All the Voigt profile fits are included in Figure 1, 
which shows the absorption-line spectra for each system in 
our sample. The total column density of each high ion and 
of H I in each system is given in Table 1. The detailed results 
of the component fits to O VI, N v, C IV, and H I in the weak 
sample are given in Tables 2 and 3. All velocities in Tables 
2 and 3 are given on a scale relative to the system redshift 

2abs ■ 

In addition to the component fits, we measured the 
total column densities in each system using the apparen t 
optical depth (AOD) method of lSavage fc Sembachl i|l99l| y 
This method has the advantage of requiring no assumptions 
about the component structure: it simply requires a choice 
of minimum and maximum velocities for the integration. It 
will return accurate column densities so long as saturated 
structure is not present. Whereas saturation is obvious if re- 
solved, because the flux goes to zero, unresolved saturation 
can be difficult to detect. In the case of a non-detection of 
Civ or Nv, we measured the (null) equivalent width in this 
ion over the velocity range in which Ovi is detected, and 
found the 3a error on this measurement. We then converted 
this maximum allowed equivalent width into a column den- 
sity limit using a linear curve-of-growth. For example, if the 
equivalent width measurement was 0±10 mA, we derived a 
3a column density limit -/Vn m using EWn m =30 mA and the 
relation N\ im = 1.13 x 10 17 EWi im /A 2 / (wi th AW i n cm" 2 
and A in A). Atomic data were taken from iMorton] (|2003l V 

In many of the strong absorbers, partial coverage of 
the continuum source is present. Indeed, we use this as a 
criterion to classify an absorber as intrinsic. Partial cover- 
age is indicated by a similar velocity structure in the two 
members of a doublet, but with inconsistent optical depths 
in each ion. In the case of partial co verage the equation 
F( v ) — Foe~ T ^ is no longer valid (e.g. iHamann fc Ferlandl 
1999), and a velocity-dependent coverage fraction must be 
introduced. For these cases, the column density derived from 
an AOD integration should be considered as a lower limit. 



3 Available at http://www.ast.cam.ac.uk/~rfc/vpfit.html 



3 RESULTS 

In this section we present a discussion of the overall proper- 
ties of the population of proximate Ovi absorbers at z=2- 
3. Although we discuss ionization processes, a detailed dis- 
cussion of the physical conditions in the individual systems 
that make up our sample is beyond the scope of this paper. 
Some of the systems in our sample are discussed in more 
detail elsewhere. These include the proximate absorber s 
toward HE 2347-4342 (|Fechner. Baade. fc Reimersl 120041 ') . 
and toward HE 1158-1 843, PKS 0329-255, and Q0453-423 
l|D'Odorico et al.ll2004h . 



3.1 Location of absorbers in velocity space 

The number of proximate O VI absorption line systems seen 
in each quasar spectrum varies between zero and five. There 
are two quasar spectra in our sample without O VI ab- 
sorption detected within 8000 kms -1 of z qBO : HE 2217- 
2818 (z qso =2.41) and Q0002-422 (z qso =2.77). There are two 
quasars with five proximate O VI absorbers in their spectra: 
PKS 0237-23 (z qso =2.233) and HE 1341-1020 (z qso =2.135). 
To illustrate the location in velocity space of the absorbers, 
we show in Figure 2 the velocity displacement of all ab- 
sorbers, strong and weak, for each quasar. Over the whole 
sample (16 spectra) we find a mean of 0.6 strong proximate 
Ovi systems per sight line, and 1.6 weak proximate Ovi 
absorbers per sight line. However, the O VI sample is incom- 
plete, because of regions where Ovi systems are missed due 
to blending by the Lya forest. As we show in §3.7, when 
we calculate the corrected redshift path for the O VI search, 
our estimated O VI completeness is 84% down to a limiting 
equivalent width of 8 mA. 

As a useful gauge of the velocity interval in which we 
expect to see proximity effects, we can derive the distance 
from each quasar at which the flux of extragalactic back- 
ground (EGB) radiation equals the flux of quasar radia- 
tion, i.e. we can derive an estimate of the size of the quasar 
sphere-of-influence. The monochromatic quasar ionizing lu- 
minosities at the Lyman Limit have been de rived for ten 
of the sixteen quasars by iRollinde et all (|2005T ). who calcu- 
lated Z/912 by extrapolating the QSO B-magnitudes assum- 
ing a spectral slope of a—— 0.5 (where F\ ~ A"). The EGB 
flux at the Lyman Limit at z=2.5 h as been calculated to be 
Fegb~6x10- 21 erg cm" 2 s" 1 Hz" 1 (|Haardt fc Madaull 19961 . 
corresponding to log Ji,=— 21.3 erg cm~ 2 s _1 Hz -1 sr _1 ; 
the values for the 2=2 and 3 cases are not significantly dif- 
ferent). The distance D at which -Fqso=Fegb is given by 
D—y^Lg^/^FEGB)- In the case where the velocities are 
Hubble flow-dominated, we can convert our calculated val- 
ues of D to velocity displacement from the quasar, using a 
value for the Hubble parameter calculated at the redshift 
of each quasar using H 2 (z) = _Hq[^m(1 + z) 3 + Q^]. With 
our adopted cosmology, //(2.5)=246 kms -1 Mpc -1 . These 
results are shown in Table 4, and are included in Figure 
2, where it can be seen that if the velocity displacements 
from the QSO are due to the Hubble flow, then the proxim- 
ity effect should only extend to ~1 000-2 500 kms" 1 . While 
peculiar velocities are to be expected, this is still a useful 
order-of-magnitude estimate for the size of the quasars' ion- 
izing sphere-of-influence. 
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Figure 2. Graphical display of the location in velocity space 
of the O VI absorbers in this study, relative to their background 
quasar. Weak absorbers are shown in red, and strong absorbers in 
blue. The location of 5^=0, i.e. z a b s =Zqso, is shown with a solid 
black line. Negative values for v qBO — v^ s imply the absorber is 
approaching the quasar. Dashed green lines show the velocity at 
which /'qso ''kgh assuming Hubble flow (see text). 



Table 4. Size of quasar sphcres-of-influence at 912 



QSO 


logLgi^ 
(erg s — 1 Hz -1 ) 


(Mpc) 


(kms -1 ) 


HE 0001-2340 


31.65 


7.7 


1700 


HE 2217-2818 


31.99 


11.4 


2600 


PKS 0237-23 


31.67 


7.8 


1700 


PKS 0329-255 


31.58 


7.1 


1800 


PKS 1448-232 


31.53 


6.7 


1400 


Q 0002-422 


31.72 


8.3 


2200 


Q 0109-3518 


31.82 


9.3 


2200 


Q 0122-380 


31.63 


7.5 


1600 


Q 0329-385 


31.28 


5.0 


1100 


Q 0453-423 


31.71 


8.2 


2100 



a Monochromatic luminosity at Lyman Limit, der ived from the 
B magnitude and F x ~ A -0 - 5 llRollinde et al.ll2005ri . 
^ Distanc e at which -Fqsq=-Fegb = 6xl0 -21 erg cm -2 s -1 Hz -1 



at z=2.5 jHaardt fc Madaul ll996l ). 



c Velocity separation from QSO at D assuming Hubble flow, with 
H(z) calculated at z qso . H(z=2.5)=246 kms -1 Mpc -1 . 



3.2 Presence of other ions 

The presence and strength of absorption in other ionic 
species in the proximate O VI absorbers provides important 
information on the ionization conditions. H I absorption is 
detected in each of the 35 O vi systems, always in Lya and 
occasionally up to Ly5. Higher-order Lyman lines may also 
be present, but they lie deep in the Lycv forest where blend- 
ing and low signal-to-noise aggravate the search for detec- 
tions. All but two of the O VI systems show C IV detections 
(the exceptions are the absorbers at z a b s =2.2378 toward 
PKS 0237-23 and z abs =2.6610 toward PKS 0329-255, both 
of which are weak O VI systems with low O VI column den- 
sities). 9/9 strong Ovi systems contain Nv, whereas 9/26 
weak O VI systems contain N v. Therefore while the presence 
of Nv is a feature of all strong systems, it does not imply 
that an O VI system is strong. Finally, C ill absorption is seen 
in 13/26 weak systems, and 7/9 strong systems. However, 
because Cm A977 lies in the Lya forest, it is often blended 
and hence not detectable even if present, so these fractions 
should be treated as lower limits. 



3.3 O VI system column densities 

A plot of the total O VI column densities in each absorber 
versus proximity (Figure 3, top panel) shows the distinction 
between the strong and weak populations. The strong pop- 
ulation all show logjV> 15.0 and are(with two exceptions) 
clustered around z a b s =2qso, whereas the weak population 
lies between the detection limit of log N ft! 12.80 (correspond- 
ing to a limiting equivalent width of 8 mA) and 14.5. Once 
the strong systems have been removed, the underlying weak 
population shows no dependence of N(0 Vl) with proxim- 
ity to the quasar. To investigate this further, we split the 
sample into Ovi absorbers at Sv < 2 000 kms -1 and those 
at 2 000 < Sv < 8 000 kms -1 . This division is motivated 
by an observed enhancement in dN/dz below 2 000 kms - 
(§3.7), as well as this representing the expected size of the 
proximity zone (§3.1). A two-sided Kolmogorov-Smirnov (K- 
S) test shows that the probability that the distribution of 
logTV(Ovi) for absorbers at <2 000 kms -1 is different from 
those at >2 000 kms -1 is only 17%. There is also no clear 
correlation between the quasar redshift and the total Ovi 
column, as shown in the lower panel of Figure 3. We do note 
that the ratio of weak systems to strong systems changes 
from 6.3 (19/3) at z < 2.5 to 1.2 (7/6) at z > 2.5, but the 
statistics are too small for us to investigate this further. 



3.4 Ovi component column densities and fr-values 

Proximity effects may occur at the component level, as well 
as the system level. To investigate this, we compare in Figure 
4 the component properties returned by VPFIT (column 
densities and line widths) of the weak Ovi absorbers with 
Sv. We include for comparison the Ovi components from 
the intervening sample of BH05. As with the system column 
densities, there is no trend for the component O VI column 
densities to correlate with proximity to the quasar. 

The distributions of log N and b for the O VI compo- 
nents in our sample are shown in Figure 5. We include the 
following component distributions as separate curves: the 
components at Sv < 2 000 kms -1 , those at 2 000 < Sv < 
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Figure 3. Scatter plots comparing the integrated column density 
(over all components) of the O VI absorbers with their proxim- 
ity to the background quasar Sv = v qBO — v^ a and with their 
redshift z a b s . An error of 500 kms" 1 in z qso is assumed. Strong 
absorbers are shown in blue (lower limits), with weak absorbers 
in red. The distinction between the weak and strong populations 
is evident. The dashed line shows the detection limit. Among the 
weak absorbers there is no correlation between JV(Ovi) and Sv, 
or between N(Ovi) and z. 

8 000 kms -1 , and the intervening absorbers from BH05. 
A two-sided K-S test finds the probability that the log N 
distributions for the samples at <2 000 kms -1 and 2 000- 
8 000 kms -1 are different (i.e. drawn from distinct parent 
populations) is only 23%. Comparing the 5f<2 000 kms -1 
components to the intervening sample, we do not find any 
difference valid at more than the l.Oer level, with the inter- 
vening sample showing a median log N(0 vi) of 13.54, and 
the Sv < 2 000 kms" 1 sample showing a median of 13.45. 
This is an important result considering the intervening sam- 
ple was measured independently in BH05, not in this anal- 
ysis. 

Turning to the component line widths, we find that 
among the proximate weak sample, 29 out of 48 components 
(60%) are narrow (6^ 14 kms -1 J3- Comparing the &-value 
distributions, we find that the probability that the popula- 
tion at Sv < 2 000 kms -1 is different from the intervening 
distribution is less than 50 per cent. The median b- value 
in the intervening sample is 12.7 kms -1 ; in the sample at 
Sv < 2000 kms -1 , it is 12.3 kms -1 . However, an excess of 
systems with 6=20-30 kms -1 is observed for components in 

4 We define 'narrow' as b < 14 kms -1 since this corresponds 
to the thermal line width of an Ovi line in gas at 188 000 K, 
and negligible Ovi is produced in gas in collisional ionization 
equil i brium (CIE) below this temperature (Sutherla nd fc Dopital 
ll993l;lGnat fc Sternber3l2007fl . 
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Figure 4. Scatter plots comparing the column densities and line 
widths of the O VI components with their proximity to the back- 
ground quasar. The errors shown are returned by the VPFIT 
profile-fitting code. Included in green is the BH05 sample of in- 
tervening Ovi components at 2=2-2.5, plotted for convenience 
at Sv > 9 000 kms -1 , with arbitrary shifts in the x-direction ap- 
plied for clarity. Two proximate components with t>=34±14 and 
45±18 kms -1 are not shown because of their large errors. 



the 2 000-8 000 kms -1 range, as can be seen in the lower 
panel of Figure 4. 

Finally, we compare the line width distribution of the 
Ovi components in the weak proximate sample with the 
C IV and N v line width distributions (in the same sample) . 
The line widths of C IV and N v distribute differently than 
those of O VI, without the strong tail extending between 10 
and 25 kms -1 seen in the Ovi distribution. 77% of Civ 
components and 90% of N v components show b < 14 km s -1 . 
These differences indicate that it cannot be assumed that the 
C IV, N v, and O VI arise in the same volumes of gas. 

In summary, among the weak O VI absorbers at z=2- 
3, there is no compelling evidence for any difference in the 
column density and line width distributions between the in- 
tervening and proximate Ovi populations, but within the 
proximate sample, the Ovi line widths distribute differently 
than those of N v and C IV. 
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Figure 5. Top two panels: normalized column density- and line 
width distributions of the proximate Ovi components at Sv < 
2 000 km s - 1 , at 2 000 < Sv < 8 000 km s ~ 1 , and in the intervening 
sample. No statistically significant differences are found, with the 
exception of a slight excess of systems with 6=20—30 kms" 1 in 
the weak sample in the 2 000—8 000 km s —1 bin. The bottom panel 
compares the line width distributions for O VI, C IV, and N V in 
the weak, proximate sample. 



3.5 H I column density 

The median and standard deviation of the logarithmic 
H I column density in the weak proximate O VI sample 
is log iV(H i)=14.66±0.95, with a total range covered in 
logiV(Hl) of over four decades. The Hl/Ovi ratio (inte- 
grated over all components in each system) varies enor- 
mously, taking values between 0.02 and 1500. In Figure 6 
(top panel) we plot the H i/O VI column density ratio against 
the Hi column density. Within both the strong and weak 
samples, a trend is evident in which Hl/Ovi rises with in- 
creasing iV(Hl). However, the two populations are clearly 
distinct, in that the weak systems show much higher H i/O VI 
ratios at a given iV(Hl). The clear segregation of the weak 
and strong systems on this plot supports our contention that 
they represent two separate populations. 

In the lower panels of Figure 6, we explore whether 
the Hi column density in the weak proximate sample de- 
pends on Sv. We find that 7V(Hl) distributes very differ- 
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Figure 6. Analysis of 7V(Hl) and the Hl/Ovi ratio in the prox- 
imate O VI sample. The column densities are integrated over all 
components in each absorber. In the top panel, a dashed line de- 
lineates the locations of the strong (blue) and weak (red) samples 
in a plot of Hl/Ovi vs Hi. In the central two panels we show 
that both ./V(Hl) and Hl/Ovi show a dependence on proximity. 
The bottom panel shows the histogram of W(Hl) for the weak 
O VI absorbers, split into two sub-samples: those at Sv < 2 000 
and Sv > 2 000 kms -1 . A significant difference exists between 
these two histograms, with the absorbers closer in velocity to the 
quasar showing lower Af(Hl). 



ently at Sv < 2 000 kms" 1 than at Sv > 2 000 kms" 1 . The 
AT (Hi) distribution for the <5u>2 000 kms" sub-sample is 
strongly clustered around log A r (H l)«15, whereas the dis- 
tribution for the Sv < 2 000 kms" 1 sub-sample is centered 
near log AT(H l)~14 (see bottom panel). A K-S test shows the 
two distributions differ at the 99.97% level. Consequently, 
there is also a difference in the mean Hl/Ovi ratio above 
and below 2 000 kms" 1 , but this is due to the proximity 
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effect in 7V(Hl), not in TV(Ovi), which shows no correla- 
tion with proximity to the quasar. The finding that iV(H i) 
correlates with Sv has an important implication: it implies 
that Sv does correlate with physical distance to the quasar. 
If Sv was dominated by peculiar motions rather than Hub- 
ble flow, and so was uncorrelated to physical proximity to 
the quasar, then the two JV(H i) distributions shown in the 
bottom panel of Figure 6 would not be distinct. The fact 
that they are so different allows us to make the general as- 
sumption that proximity in velocity does imply proximity in 
distance. 

3.6 Civ and Nv column densities 

The logarithmic C IV column densities in the proximate O VI 
sample cover over 2 dex, with a median value and standard 
deviation of log N(C iv) = 12.64±0.60. N v is detected in nine 
weak absorbers, with logiV(Nv) between 11.54 and 13.60. 
Comparisons between the Civ column density and Sv, and 
between Nv column density and Sv, as well as the depen- 
dence of the high-ion column density ratios with Sv, are 
presented in Figure 7. In a similar fashion to Hi, iV(Clv) 
in the weak Ovi absorbers depends strongly on Sv, with 
the absorbers at Sv < 2 000 kms" 1 showing column densi- 
ties clustered around log iV(C iv)~12.4, and those at Sv > 
2 000 kms -1 showing values centered at log iV(C iv)~13.3. 
A two-sided K-S test shows that the iV(Clv) distributions 
in the two velocity intervals differ at a significance level of 
93%. This leads to a difference in the mean Civ/Ovi ratio 
between the two velocity bins, with lower ratios at lower ve- 
locity. There is no indication that either jV(N v) or N v/O VI 
depends on Sv, though the smaller sample size of Nv (be- 
cause of the non-detections) prevents us from making sta- 
tistically meaningful conclusions. 

In order to study the high-ion velocity structure more 
closely, we present in Figures 8a and 8b apparent column 
density profiles as a function of velocity, for each weak 
proximate O VI absorber. The apparent column density per 
unit velocity in units of ions cm " 2 (kms" 1 )" 1 is given b y 
N a {v) = 3.768 x 10 14 T a (v)//A |s avagc fe Sembach|[l99l l. 
where A is the transition res t wavelength in A, / is the os- 
cillator strength (taken from lMortonll200"^ 'l . and T a (v) is de- 
fined in §2.2. These figures allow the Ovi, Civ, Hi, and 
(in the nine cases where present) Nv profiles to be closely 
compared. Note that for Hi, we plot the apparent column 
density profile of our best-fit model to the Lyman series ab- 
sorption lines, because several lines are used to derive our 
best-fit model. For the high ions, the actual data from one 
of the two doublet lines are shown, rather than the model. 

In 11/24 absorbers where both Ovi and Civ are 
detected, the profiles are similar enough to suggest co- 
spatiality of the ions C +3 andO +5 .In the remaining 13 cases 
there are significant differences in the line profiles, indicating 
that Ovi and Civ trace different gas phases. Here a 'signif- 
icant difference' is either a centroid offset of >10 kms"" 1 , or 
a case where &(Civ)<6(Ovi), for which there is no single- 
phase solution [since an oxygen atom is heavier than a car- 
bon atom, &(Ovi) must be less than 6(Civ) for lines arising 
in the same gas] . In only one of the nine cases with N v, its 
profile is similar in centroid and width to the O VI. How- 
ever, detailed comparisons of the N v and O VI profiles are 
complicated by the fact that low-optical depth Nv absorp- 



Table 5. Proximate Civ with blended Ovi 



QSO 






6v (kms x ) 


HE 0001-2340 


2.2670 


2.1870 


7530 


HE 0151-4326 


2.7890 


2.6956 


7580 


HE 1341-1020 


2.1350 


2.1065 


2750 


HE 2347-4342 


2.8710 


2.8781 


-550 


Q 0122-380 


2.2030 


2.1472 


5320 



tion is lost in the noise, even in this high-quality dataset, 
so we do not read into this result any further. Most sig- 
nificantly, offsets of order 10-30 kms - between the cen- 
troid of the H I absorption and the centroid of O VI absorp- 
tion are seen in 14/26 weak proximate cases. These off- 
sets are highly important, since they imply that the H° 
atoms do not live in the same volume of gas as the +5 
ions in at least half of the weak O VI absorbers. In other 
words, at least half of the weak O VI systems are multi- 
phase. Consequently, it is not possible to derive a metal- 
licity estimate in these multi-phase cases from the Ovi/Hl 
ratio, even with an ionization correction, and single-phase 
photoionization models cannot be applied. The dangers of 
treating multi-phase qu asar absorbers as single-phase have 
been pointed out be fore l|Giroux. Sutherland, fe Shulllll994l ; 
iReimers et ai]|200ll ). 

3.7 Calculation of dN/dz for proximate Ovi 

In this section we derive the incidence of proximate O VI 
absorption, in terms of the traditional dN/dz, and also in 
terms of the number of absorbers detected per 2 000 kms" 1 
interval in Sv. Note that because we are working at z 3> 0, rel- 
ativistic effects cause the relationship between velocity and 
redshift intervals to be dv « cdz/(l+z). Before calculating 
the incidence of proximate Ovi absorption, we correct the 
redshift path for blending with the Lya forest, since O VI is 
only identifiable in unblended regions of the spectrum. To 
estimate the number of O VI absorbers that were missed due 
to Lya forest contamination, we looked for C IV absorbers 
within 8 000 kms" 1 of z qBO in which the accompanying Ovi 
is blended, and found five. These are listed in Table 5. This 
method should provide a good estimate for the O VI incom- 
pleteness, since firstly O VI is detected in all proximate C IV 
absorbers in which the O VI data are unblended, and sec- 
ondly C IV is unblended by the Lya forest, so no C IV sys- 
tems should be missed down to a limiting equivalent width 
of «2 mA. The limiting equivalent width for the C IV search 
is smaller than the limiting equivalent width for the O VI 
search, because C IV lies in a higher S/N region of the spec- 
trum. 

The total redshift path within 8 000 kms" 1 of sixteen 
quasars (at (z qso )=2.44) is 1.48. Because there are five 'miss- 
ing' Ovi systems, we estimate our Ovi completeness to 
be 35/ (35+5) =84%, implying that a path equal to 0.25 is 
blended near O VI, so that the corrected redshift interval for 
the proximate Ovi search is 1.24. In this path 26 weak sys- 
tems (comprising 48 components) are detected, so overall 
dN/dz of weak Ovi systems at Ji><8 000 kms" 1 is 21±5, 
where the errors on the counts are assumed to be Poissonian. 
However, this is an average value, and when we break down 
the absorbers into those at Sv < 2 000 and those at 2 000- 
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Figure 7. Dependence of ./V(0 vi), N(C iv), 7V(C rv)/iV(0 vi), N(N v), and JV(N v)/N(0 vi) on Sv, for the weak proximate O VI sample. 
Downward-pointing arrows denote non-detections. The right column shows the histograms of each quantity plotted on the left, split 
into sub-samples at Sv < 2 000 and >2 000 kms" 1 from the quasar. In the lower two panels, the nine Nv detections are shown with 
open circles to distinguish them from the non-detections. Only detections are included in forming the histograms. Whereas there is no 
discernable effect for ./V(Ovi) to depend on Sv, JV(C iv) distributes very differently at <2 000 than at >2 000 kms -1 . The Nv samples 
arc too small to draw conclusions about whether N(Nv) depends on Sv. 



8 000 kms , we find a marked change in the incidence: 
diV/dz trebles within 2 000 kms -1 , from 14±4 to 42±ia3. 
When counting components rather than systems, a similar 

5 We have included the four weak O VI absorbers at Sv < 
(z a bs>Zqso) i n the 0—2000 kms - bin, since regardless of their 
peculiar velocity they must lie in front of the quasar. Even with- 



increase in the incidence is observed within 2 000 kms -1 of 
z qso . These results are shown in Figure 9, and are also sum- 
marized in Table 6. Note how the intervening incidence of 
BH05 is recovered in the range 2 000-8 000 kms -1 . 

out these, dN/dz in the 0-2 000km s -1 bin is 28±8, twice the 
incidence in the 2 000-8 000 kms -1 bin. 
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Figure 8. High-ion and Hi apparent column density profiles for each weak proximate Ovi absorber. The Civ (blue), Hi (black), and 
Nv (green, where present) profiles have been scaled to allow their shape to be compared with Ovi. For the high ions the AOD profiles 
are formed from the data themselves. For H I, we show the AOD profile of the best-fit model to the Lyman series absorption lines. The 
green vertical lines show the velocity range over which the AOD measurement is made. Annotated on the top-right of each panel is the 
proximity to the quasar. In cases where the Ovi and Civ profiles show similar line centers and line widths, we add a tick to the Civ 
label. In cases where the Ovi and Hi profiles show similar velocity centroids (but not necessarily the same width), we add a tick mark 
to the Hi label. Cross marks next to the Civ and Hi labels indicate significant differences between the profiles of these ions and Ovi, 
e.g. centroid offsets of >10 kms -1 , indicating that the lines do not form in the same phase of gas as the Ovi. 



3.8 Comparison of z ~ 2.5 and z ~ absorbers 

A considerable amount of work has been invested in studying 
O vi in the low-redsh ift IGM using the FUSE and HST or- 
bitin g observatorie s (iBergeron et al.lll994l; iBurles fc Tvtlerl 
Il996l ; iTripp et all Il996l ; iTripp. Savage, fc Jenkins! I2OO0F 



Savage et al. 20021; ISembach et al. 20041; Richter et al. 


2004; Prochaska et all 12004; 


Danforth & Shulll 12005, 


2008; 


Lehner et al. 2006; Tripp et all 2008; Cooksev et all 


2008; 


Thorn & Chen 2008a|lJ. By comparing our results obtained 



from proximate O vi absorbers at z=2-3 with the 2 w re- 
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Figure 8. (continued). High-ion and Hi apparent column density profiles in the weak Ovi absorbers. 



Table 6. Weak Ovi absorbers at 2=2-3: Summary of dN/dz 
statistics (8 mA sensitivity) 



Category | <5u | < 2 OOO^ 2 000 < Sv < 8 000 a Intervening^ 

Systems 42±12 14±4 14±2 

Components 87±17 23±5 35±6 



a Range in velocity relative to z qBO (in km s 1 ) . 
" Based on BH05; in the updated intervening sample 58 sys- 
tems are seen over a total Az of 4.09, complete down to 
logjV(Ovi)=13.11. 

suits, we can search for evolution of the highly ionized gas 
content of the Universe. Since the intensity of the extragalac- 
tic U V background rises with redshift (H aardt fc Madaul 
1996), one might expect to see higher photoionized O vi col- 
umn densities at higher redshift. Simultaneously, the frac- 
tion of baryons in w arm-hot gas is pre dicted to decrease with 
increasing redshift ((5 ave et al.ll200lD . since at early cosmic 
epochs there has not been adequate time for gas to accrete 
and virialize into dark matter-dominated potential wells. 



Based on an analysis of low-redshift (0.15 < z < 0.5) 
Ovi absorption line systems and w ith a sensitivity limit 
of W r > 30 mA, iTripp et all (|2008l ) report that dN/dz is 
enhanced by a factor of ~3 at velocities <2 500 kms" 1 
from the quasar compared to the intervening sample, 
with dN(0 Vl)/di?=16!2 f° r the intervening systems and 
diV(Ovi)/d2=5lli5 for those at Sv < 2500 kms" 1 . Our 
result for weak O vi systems at z=2-3 with W T > 30 mA 
and Sv < 2 000 kms" 1 is dN/dz is 32±10, lower than the 
low-redshift result by a factor of ~1.6. However, comparing 
the high- and low-redshift results for dN/dz can be mis- 
leading, because redshift has a highly non-linear relation- 
ship with spatial dimensions, that depends on cosmology. 
With our adopted cosmology, then in the interval 2=0.15 
to 2=0.5, d:r(Gpc)=3.63 dz, whereas for the interval z=2 
to z=3, dx(Gpc) = 1.22 dz, where da: is the comoving radial 
distance, and where we have used the online cosmology cal- 
culator of lWrighd 1120061 ) . Correcting for these factors to cal- 
culate the number density of weak O vi absorbers per unit 
co-moving distance, a better indicator of the space density 
of Ovi, we derive the results listed in Table 7. Here we see 
evidence that the space density of Ovi absorbers increases 
with redshift, by a factor of two between z=0.15-0.5 and 
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Figure 9. Number of weak Ovi absorbers as a function of prox- 
imity to the quasar. Plotted are the number of O VI systems per 
2 000 km s — 1 bin of Sv over our sample of sixteen quasars. The cor- 
responding values of dN /dz are shown on the right, corrected for 
the blended redshift path. The errors on the counts are assumed 
to be Poissonian; plotted are the la errors. We use a smaller bin 
for the range —1000 to kms -1 ; these absorbers lie in front 
of the quasar and so are included in the 0-2 000 kms -1 bin in 
our final dN/dz calculations. The incidence of intervening Ovi 
systems is shown as the green data point; there is no significant 
difference in dN/dz between the absorbers at 2 000—8 000 kms 
and the intervening sample. 



Table 7. dA r (Ovi)/dx statistics (30 mA sensitivity) 



Category^ 


Sample 


z=0-0.5 b 
(Gpc- ~ 


2=2-3 
(Gpc" 1 ) 


Systems 


Proximate (Sv < 2 000) 
Intervening (Sv > 5 000) 


14±« 
4±1 


26±8 
8±2^ 


Components 


Proximate (Sv < 2 000) 
Intervening (Sv > 5 000) 


30t 8 6 
6±1 


71±14 
21±3 C 



a The entries in this table list the incidence of weak Ovi ab- 
sorbers per unit comoving distance, in units of Gpc -1 ; the en- 
tries were derived by correcting dN/dz for cosmology, and are 
presented for both proximate and intervening absorbers, at low 
and high redshift, and for systems and components. 
" Derived using results from lTripp et all |2008). 
c Derived from BH05 result. 

2=2-3. This increase is observed independently for both the 
proximate and intervening absorbers. 



4 DISCUSSION 

The BAL, mini-BAL, and intrinsic systems, which together 
we classify as strong systems based on their fully saturated 
O VI absorption, strong accompanying O VI and C IV absorp- 
tion, and frequent evidence for partial coverage, are well- 
understood as being formed in either QSO-driven outflows 
or inflows in the i mmediate vicinity of the A GN central en- 
gine (see review bv lHamann fc Ferl and 1999). Here we focus 
on the weak (narrow) proximate systems, which represent a 



separate population. Is there any evidence that these weak 
systems are directly photoionized by the quasar? 

We report two results that address this question, though 
their interpretation is not straightforward. First, in our sam- 
ple of proximate O VI absorbers at 2=2-3, there is an en- 
hancement by a factor of three in the incidence dN/dz of 
weak Ovi absorbers within 2 000 kms -1 of the quasar ver- 
sus those in the interval 2 000 to 8 000 kms - 1 . A similar 
enha ncement is seen near low-redshift quasars ([Tripp et al.l 
2008). While at face value this could be interpreted as 
a proximity effect, in which quasars preferentially ionize 
nearby clouds more often than they photoionize more dis- 
tant clouds, this is not the only explanation. The en- 
hancement in dN/dz could also be explained by an over- 
density of galaxies near quasars, and where the Ovi ab- 
sorbers are located in the gaseous h alos of these galaxies 
(| Young. Sargent, fc BoksenbergHl982r ) . It is well-known that 
quasars are preferentially fo rmed in cluster environments, 
parti cularly at high redshift |Croom et alj |2005; She n et all 
l2007h . This idea is supported by the observation that the 
internal properties of the weak O VI absorbers are not de- 
pendent on redshift or proximity to the quasar. 

Second, there are statistically significant differences be- 
tween the Hi/Ovi and Civ/Ovi ratios measured in the 
weak Ovi populations above and below 2 000 kms -1 , with 
the weak absorbers at Sv < 2 000 kms -1 showing a me- 
dian Hi/Ovi ratio lower by a factor of ~10, and a me- 
dian Civ/Ovi ratio lower by a factor of «7 than those at 
2 000 < Sv < 8 000 kms -1 . A gain, at face value, these results 
could be interpreted as proximity effects, in which gas closer 
to the quasar shows a higher ionization level than interven- 
ing gas. However, closer inspection finds that while iV(Hi) 
and AT(Civ) show a tendency to decrease as Sv decreases, 
N(Ovi) is uncorrelated with Sv, i.e. it is the behaviour of 
iV(Hi) and AT(Civ) alone that is driving the trends seen 
in the column density ratios, not the behaviour of JV(0 Vi). 
Importantly, the H I and C IV proximity effects imply that 
proximity in velocity does correlate with proximity in dis- 
tance, which supports the idea that the velocities of the weak 
absorbers are dominated by the Hubble flow. 

Because Ovi-Hi-Civ velocity centroid offsets are ob- 
served directly in «50 per cent of the weak systems, single- 
phase photoionization models for the Ovi, Civ, and Hi 
are inadequate for at least half the absorbers in our sam- 
ple. Indeed, the non-dependence of N(0 vi) on proxim- 
ity casts doubt on whether photoionization by the back- 
ground QSO creates the Ovi at all. Photoionization by 
nearby stellar sources of radiation is also unlikely, since 
such sources do not emit sufficient fluxes of photons 
above 54 eV (the Hen ionization edge) to produce the ob- 
served quantities of O VI. We cannot rule out photoion- 
ization by the quasar in every individual case, for exam- 
ple the absorber at 2=2.4183 toward HE 1122-1649, which 
shows strong Ovi but barely detectable Civ, as do sev- 
eral p roximate O vi absorbers reported by iGoncalves et al.l 
l|2008l ). Nonetheless, JV(Ovi) does not depend on proxim- 
ity in the way that A r (Hl) and iV(Clv) do. Thus we in- 
fer that the sizes of spheres-of-infiuence around quasars, 
which are implied the shapes of Gunn-Peter sen troughs 
([Zheng fc Davidsenl 1 19951 ; ISmette et all |2002|) . and from 
studies of proxi mate absorption in lower ionization species , 
such as Mg ii (|Vanden Berk et all 120081 ; IWild et alj l2008h 
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and Civ l|Foltz et al.lll986l ; IVestergaardll2003l ; I Nestor et al.l 
I2008T I. are dependent on photon energy, and the presence of 
a sphere of photons at energies above 113.9 eV (capable of 
ionizing O" 1 " 4 to +5 ) has yet to be demonstrated. 

Thus we turn to collisional ionization models. Since 
the line widths of a significant fraction (~60%) of the 
components in the weak O VI absorbers are low enough 
(fo<14 kms -1 ) to imply gas temperatures below 188 000 K, 
collisional ionization equilibrium can be ruled out, because 
essentially no O VI is produced in gas in CIE at these tem- 
perat ures (jSutherland fc Dopital 1 19931 ; iGnat fc Sternberg! 
120071 ). Indeed, the narrow line widths of many intergalac- 
tic Ovi absorbers at z ~ 2 have led various authors 
to conclude that photoionization is the origin mechanism 
llCarswell et~a l. 2002; B ergeron et~alir2 002: iLevshakov et ajj 
20031: Bergeron fc Hcrbcrt-Fort 2005; Rci mers et all l200rj ; 



Lopez et al.ll2007f ). However, non- equilibrium collisional ion- 
ization models cannot be ruled out so easily. Indeed one 
expects that collisionally ionized gas at 'coronal' tempera- 
tures of a few x 10 s K, where O VI is formed through colli- 
sions, will be in a non-equilibrium state. This is because the 
peak of the interstellar cooling curve exists at these temper- 
atures, and so the cooling timescales are short. When the 
cooling times are shorter than the recombination timescales, 
'frozen-in' ionization can result at temperatu res well below 
those at which the ions exist in equilibrium dKafatosI 1 19731 ; 
IShapiro fc Moore! Il976l ; lEdgar fc Chevalier! Il986l ). provided 
that there is a source of ~ 10 6 K gas in the first place. 

There are at least two physical reasons why 
collisionally-ionized, million- degree regions of interstellar 
and intergalactic gas could arise in the high-z Universe. 
The first is the (hot-mode) accretion and shock heat- 
ing o f gas falling into p otential wells (jBir nboim fc Dekell 
120031 ; iKeres et al] 120051 ; iDekel fc Birnboiml l200fj ). a pro- 
cess which is inc orporated into cosmological hydrodynam- 
ical simulations llCen fc Ostrikei 1 19991; iDave et all [200H ; 
iFang fc Brvanll2003l ; iKang et al.l [20051 : ICen fc Fang||2006h . 
and creates what is referred to as the Warm-Hot Inter- 
galactic Medium. However, these models generically predict 
that the fraction of all baryons that exist in the WHIM 
rises from essentially zero at z=3 to 30-50% at z=0, so 
little WHIM is expected at the redshifts under study here. 
The second reason is the presence of galactic-scale out- 
flows, which due to the energy input from supernovae are 
likely to contain (or even be domi nated by) hot, highly 
ioniz e d gas (see recent models by lOppenheimer fc Dave 



2006:; iFangano. Ferrara. fc Richterl 1200 
2007; ISamui. Subramanian. fc Srianan 



002 

and 



Kawata fc Rauch 
2008). There 



strong observational ev idence for outflows at redshifts of 
«2-3 |Heckmanl I2002T I. including blueshifte d absorption 
in the spe c tra of Lyman brea k galaxies (|Pettini et al.l 
l200d . 120021 ; IShaplev et ail 120031 ). the presence of metals 
in the low-density , photoionized I GM ( the Ly q forest; e.g. 
Aracil et al.l 12004 lAguirre et all [20051 . I2008T I. though see 



Schave. Carswell. fc Kiml i|2007l ). and the presence of super- 



escape velocity C IV c omponents in the spectra of damped 
Lya (DLA) galaxies (|Fox et al.|[2007br i. In addition, colli- 
sionally ionized gas in galactic halos at high redshift has 
been seen directly th rough detection s of O VI and N v 
components in DLAs l|Fox et al.l l2007ah . albeit with much 
broader system velocity widths than in the proximate ab- 
sorbers discussed here. O VI absorbers with lower H I col- 




-2 -1 
log N(N V)/N(0 VI) 

Figure 10. Comparison of predictions of non-equilibrium colli- 
sional ionization models with the observed column density ra- 
tios in our weak proximate sample. Non-detections of N V (C iv) 
lead to upper limits on the x- (y-) axis. We take the predictions 
from the isobaric non-eq uilibrium collisional ionization model of 
IGnat fc Sternberg! [120071) for an absorber at T=10 5 K (higher 
temperatures are ruled out by the component line widths). The 
model is fairly insensitive to changes in the overall abundance 
level, but is sensitive to non-solar relative abundances. The dashed 
green lines show the range of values of [N/O] and [C/O] that are 
required to fit the proximate data with this collisional ionization 
model. 



umn density may probe the outer reaches of such halos or 
'feedback zones' (BH05: [Simcoe et alJliool ). 

We explore the ability of non-equilibrium collisional 
ionization models to reproduce the data in our proxi- 
mate Ovi sample in Figure 10, which shows the Civ/Ovi 
vs Nv/Ovi ratio-ratio plane. We take the isobaric non- 
equilibrium model at logT=5.00 (consi stent with the ob- 
served Ovi component line widths) from lGnat fc Sternberg! 
(2007), computed using solar abundances, and then find the 
values of [N/O] and [C/H] that are required to reproduce 
the observations. The models for 0.1 and 0.01 solar abso- 
lute abundances, and for the isochoric case, predict similar- 
values. The non-equilibrium models are not capable of re- 
producing the observed ratios when solar relative elemental 
abundances are used. However, if [N/O] takes values be- 
tween — 1.8 and 0.4, and [C/O] between —1.9 and 0.6, it is 
possible to explain the observed column densities in each 
weak proximate system with a non-equilibrium collisional 
ionization model. Note that we have not corrected these 
model predictions for the effect of photoionization by the 
extragalactic background; the model represents the pure col- 
lisional ionization case. Hybrid collisional+photo-ionization 
models w ould help in makin g progress in this area (see dis- 
cussion in lTripp et al.|[2008l '). 



5 SUMMARY 

We have presented a study of proximate (Sv < 8 000 kms -1 ) 
Ovi absorbers in the spectra of sixteen quasars observed 
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at high signal-to- noise and 6.6 kms resolution with 
VLT/UVES. The quasars are at redshifts between 2.14 and 
2.87. We found 35 proximate O VI absorbers comprising 101 
components. We used component fitting software to deter- 
mine the properties of the absorption, studied the statistical 
properties of the sample, and addressed the ionization con- 
ditions in the gas. Our study has produced the following 
results. 

(i) Nine of the 35 Ovi systems are strong, with 
logJV(Ovi) > 15.0. These absorbers all show detections of 
Civ and Nv, and show either broad, fully saturated, mini- 
BAL-type Ovi absorption troughs extending over tens to 
hundreds of kms , or evidence for partial coverage in the 
high-ion profiles. A mean of 6.0 O VI components are seen 
in the strong absorbers. The strong absorbers are formed 
either in QSO-driven outflows or in the immediate vicinity 
of the quasar. 

(ii) 26 of the 35 O VI systems are weak, and show 
log N (Ovi) < 14.5. These systems are narrow, with a me- 
dian velocity width Avgo of 42 kms -1 . The weak systems 
all contain detectable Hi, all but two show Civ, and 9/26 
show N v. We find a mean of 1.9 O VI components per weak 
system. Approximately 60 per cent of these O VI compo- 
nents have Doppler b- values <14 kms -1 . The gas in these 
narrow components is constrained by the line widths to be 
at T< 188 000 K. 

(iii) Among the weak sample there is no correlation be- 
tween the O VI column density and the proximity to the 
quasar. This is true when using either the component-level 
or the system-level column densities. Two-sided K-S tests 
confirm that the column density distributions of the inter- 
vening and weak proximate O vi samples are statistically 
indistinguishable. The O VI column density also does not 
correlate with redshift. 

(iv) There is a difference between the line width distribu- 
tions of the O vi components versus the distributions of C IV 
and N V, in the weak proximate sample. The C IV and N V 
distributions are strongly peaked below 6=10 kms -1 and do 
not show the tail of absorbers with b- values between 10 and 
25 kms -1 seen in the Ovi distribution. 

(v) There is a significantly lower 7V(H i) in the weak 
Ovi absorbers with Sv < 2 000 kms -1 than in those with 
<5«>2 000 kms -1 . The difference in the median iV(Hl) be- 
tween these two samples is ~1.0 dex. We interpret this as 
a manifestation of the Hi proximity effect, and it leads to 
a trend in which the Hl/Ovi ratio is lower for absorbers 
within 2 000 kms -1 of the QSO than for those at higher 
Sv. An identical trend is observed for C IV, with the median 
TV(CIV) lower by ^0.8 dex at Sv < 2 000 kms -1 than at 
Sv > 2 000 km s~ 1 . Thus C IV and H I behave differently than 
Ovi. 

(vi) Apparent column density profile comparisons show 
that O vi and C IV show significantly different velocity struc- 
ture in 13 of the 24 weak O vi absorbers where C iv is de- 
tected. Offsets between the centroid of the H I model and the 
Ovi components are seen in 14/26 weak proximate cases. 
Therefore at least half of the weak Ovi systems are multi- 
phase. 

(vii) Down to a limiting rest equivalent width of 8 mA, 
the number density of weak O vi systems within 2 000 km s _1 
of the background quasar is diV/dz=42±12 (87±17 for com- 



ponents). Between 2 000 and 8 000 kms -1 in Sv, dN/dz falls 
to 14±4 (23±5 for components), equal to the incidence of 
intervening O vi measured by BH05. Thus an enhancement 
in dN/dz by a factor of three is seen at <5u<2 000 kms -1 . 

(viii) We compare the incidence dN/dz of Ovi absorbers 
at high (2=2-3) and low (2=0-0.5) redshift. After cor- 
recting for cosmology, the incidence of O vi systems with 
W T > 30 mA per unit comoving distance at z=2-3 is twice 
the value at 2=0-0.5. This increase is seen for both proxi- 
mate and intervening absorbers. 

(ix) Because O vi-C iv-H I velocity offsets are observed in 
half the sample, single-phase photoionization models are (at 
least for these cases) unable to explain the data. Indeed, 
the non-dependence of N(0 vi) on proximity casts doubt 
on whether photoionization creates the O vi at all. Instead, 
we propose that the weak O vi absorbers trace collisionally- 
ionized regions that exist as part of multi-phase galactic or 
protogalactic structures. In order to reconcile the collisional 
ionization hypothesis with the narrow line widths measured 
in approximately 60 per cent of the O vi components, the gas 
must be out of CIE. This can be explained by a scenario in 
which initially hot gas is cooling down through the coronal 
regime at T~10 J K. The observed high-ion column density 
ratios and line widths in the proximate O vi absorbers can be 
explained by models of non-equilibrium collisionally ionized 
gas at T w 10 5 K, but only if the relative elemental abun- 
dances are non-solar: we require [N/O] to be in the range 
— 1.8 to 0.4 and [C/O] between —1.9 and 0.6, to explain the 
data with these models. 

(x) Future O vi surveys do not need to remove absorbers 
(and thus compromise the sample size) that happen to lie at 
2 000-8 000 kms -1 relative to the quasar. A better approach 
would be to remove those showing either iV(0 vi)>15.0 or 
evidence for partial coverage. 
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Table 2. Component fit results to weak Ovi absorbers 






Table 3. Component fit results to weak O VI absorbers (cont.) 
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Table 1. Proximate Ovi absorbers in the UVES Large Programme 
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a Velocity separation between absorber and quasar, in kms 1 . A negative value indicates the absorber is 
at higher redshift than the quasar. Errors are ~500 kms -1 , reflecting the uncertainties in Zq ao . 
" Ovi line width containing the central 90 per cent of the integrated optical depth, in kms -1 . Only given 
for the weak absorbers, since it cannot be measured for a saturated line. 

c The column densities in this table (with N in cm -2 ) were determined by summing the column densities 
of the individual components in our Voigt profile fits, with the component errors added in quadrature. 
For saturated absorbers, we present lower limits on the column density. The symbol | denotes absorbers 
with evidence for partial coverage of the continuum source; these column densities should be considered 
lower limits. For non-detections of C IV or N V, we present a 3<r upper limit to log N, based on the null 
measurement of the equivalent width. Absorbers marked f are not formally detected at 3cr significance, but 
a fit was executed nonetheless because of weak features in the line profiles at exactly the same redshift as 
O VI (see Fig 1). 



High-z Proximate Ovi 19 



STRONG, 5v= 7227 km s J 
HE 0151-4326 z =2.< 




^ 

P 

b 1 



,0 
u 

< 



£ 4 



STRONG, <5v= 6103 km s"' 
HE 0151-4326 z =2.7134 



_ g _ VI 103 



STRONG, <5v= 557 km s~' 
HE 1158-1843 z =2.4426 



STRONG, <5v= 1740 km s~' 
HE 1311-1020 z =2.1169 




3 - 



4 - 





£.6 



D.O 

4 









H I 972 (Lyy) 




1025 (Ly|3) 



1215 (Lya) 





E 83 



H I 972 (Lyy) 



4 



H I 1025 (Ly/S) 





n u i i I i i i I i^yj I i i i I i i i I i i U o U i i I i i i I i i i I N-t-U-< q u i i I i i i I i i Tt i i i I i i i I i i u o 

-200-100 100 800 -200-100 100 200 -200-100 100 200 -800 -400 400 800 

v (km s -1 ) v (km s -1 ) v (km s" 1 ) v (km s" 1 ) 



Figure 1. Absorption-line spectra of all proximate Ovi systems in the UVES Large Programme. Strong absorbers are shown first, 
followed by the weak absorbers. Our Voigt profile fits are shown in red. In many of the C IV and N V panels, the data are shown over a 
limited range in the y-axis, for clarity. Green vertical lines indicate the velocity range of the AOD integration. The label 'PC indicates 



20 Fox, Bergeron, & Petitjean 



STRONG, <5v=-1171 km s H 
HE 1341-1020 z =2.1473 



STRONG, <5v=-1586 km s~ 
HE 2347-4342 z =2.B916 



STRONG, dv=-2015 km s H 
HE 2347-4342 z =2.8972 



STRONG, <5v = -2541 km s~ 
HE 2347-4342 z =2.9041 
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High-z Proximate Ovi 21 



STRONG, <5v= 2190 km s" 1 
PKS 0320-255 z =2.70B9 



WEAK, <5v= 911 km s" 
HE 0001-2340 z =2.2571 



WEAK, dv= 616 km s"' 
HE 0001-2340 z =2.2603 



WEAK, <5v= 5205 km 
HE 1122-1649 z =2.351B 




P 2 



T 



_ C III 977 



H t 978 (Ly 7 ) 





O VI 10: 



N V 1238 



o vi 10; 




N V 1238 



10 




VI 1031 1 1 



=H 



972 <Ly 7 ) 




1215 (Lya) 





t 



H 1 972 (Ly 7 ) 



1025 (Ly/S) 





n u i i I i i i i i m l-l / lllllljjj a Li i i i i i i i i uj^i i I i i i I i i y o yjjjjjjjjjjjjjjjjjjjjjjj o 

-200-100 100 200 -200-100 100 200 -200-100 100 200 -200-100 100 200 

v (km s"') v (km s"') v (km s"') v (km s"') 



Fox, Bergeron, & Petitjean 



WEAK, <5v= 3720 km s" 1 
HE 1122-1649 z =2.3682 



WEAK, <5v= -727 km s~ 
HE 1122-1649 z =2.4183 



WEAK, dv= 1845 km s" 1 
HE 1158-1843 z =2.4279 



WEAK, <5v= -451 km s~ 
HE 1158-1843 z =2.4542 
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High-z Proximate Ovi 23 



WEAK, <5v= 4858 km s" 1 
HE 1341-1020 z =2.0850 



WEAK, <5v= 3213 km sT' 
HE 1347-2457 z =2.5727 



WEAK, 5v= 3069 km s H 
HE 1347-2457 z =2.5744 



WEAK, <5v= 659 km s~' 
HE 2347-4342 z =2.8625 
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Fox, Bergeron, & Petitjean 



WEAK, <5v= 2826 km s 
PKS 0237-23 z =2.2023 



WEAK, <5v= 1819 km s" 
PKS 0237-23 z =2.2135 
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High-z Proximate Ovi 25 



WEAK, <5v= -444 km s" : 
PKS 0237-23 z =2.2378 



WEAK, <5v= 7114 km s~ 
PKS 0329-255 z =2.6494 
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Fox, Bergeron, & Petitjean 



WEAK, <5v= 246 km s H 
Q 0109-3518 z =2.4012 



WEAK, Sv= 4204 km s 
0122-380 z =2.1587 



WEAK, dv= 7870 km s" 1 
Q 0329-385 z =2.3520 
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High-z Proximate Ovi 27 



WEAK, <5v= 5882 km s" 1 
Q 0339-385 z =2.3738 



WEAK, <5v= 1797 km s~ 
Q 0453-423 z =2.6362 



WEAK, 5v= 1465 km s" 1 
Q 0453-423 z =2.6402 
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